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In the totally asymmetric simple exclusion process (TASEP) two pro- 
cesses arise in the large time limit: the Airyi and Airy2 processes. The 
Airy2 process is an universal limit process occurring also in other models: 
in a stochastic growth model on 1 + 1-dimensions, 2d last passage percola- 
tion, equilibrium crystals, and in random matrix diffusion. The Airyi and 
Airy2 processes are defined and discussed in the context of the TASEP. 
We also explain a geometric representation of the TASEP from which 
the connection to growth models and directed last passage percolation is 
immediate. 



(N 
O 



Abstract 



1 Introduction 
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^p H , In mathematics and physics, there are many situations where an observable 

' behaves, in an appropriate limit, as a Gaussian random variable. The most 

famous example in mathematics is arguably the central limit theorem for the 
sum of i.i.d. random variables provided the existence of more than two moments. 
This is the simplest example of universality class and extends in some cases to 
m ^ ' dependent random variables too. 

, The Gaussian is however not the only universality class. In the last decade 

;_h ' a lot of progress in understanding universality related to random matrix en- 

sembles have been achieved. In 1994 Tracy and Widom determined the limit 
distribution of the largest eigenvalue in the Gaussian Unitary Ensemble (GUE) 
of random matrices [50], thus called GUE Tracy- Widom distribution, Then, 
in 1999 Baik, Deift, and Johansson proved that the same distribution asymp- 
totically describes the longest increasing subsequence of a random permuta- 
tion [3]. Johansson figured out that F 2 occurs also in the shape fluctuation 
in a point-to-point directed last passage percolation model [28]. The same 
limit distribution arises in a stochastic growth model on a one-dimensional sub- 
strate, the polynuclear growth model [40], in other corner growth models [23], 
directed percolation [6,8,28], in vicious random walks [36], in non-colliding 
Brownian motions [32] , and in the totally asymmetric simple exclusion process 
(TASEP) [4,28,37]. Concerning random matrices, it was meanwhile proven that 
F 2 occurs beyond the GUE ensemble [13,47]. 
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In several of the models analyzed, a natural question in an extension to 
joint distributions. In random matrices, the extension is Dyson's Brownian 
Motion [14]. The process describing the evolution of the largest eigenvalue in 
(3 = 2 Dyson's Brownian Motion is, in the limit of large matrix dimension, 
the Airy process (below referred as Airy2 process). This process was discov- 
ered by Prahofer and Spohn in a growth model of a surface, the polynuclear 
growth (PNG) model [42], where the Airy 2 process describes the interface for 
large growth time t. As expected, it arises also in discretized versions of the 
model [25,26,29]. The PNG surface can be also related to an evolution of certain 
Young tableaux, in terms of which the result has been extended by Borodin and 
Olshanski [11]. Discrete versions of the PNG model are tightly linked with a 
point-to-point directed polymer problem (or directed last passage percolation) 
where again the Airy 2 process is a limit process [29]. A different class of mod- 
els where the Airy 2 process arises are tiling models, the most studied one is 
the Aztec diamond [30], which can be mapped to the six- vertex model at the 
free-fermion point [21]. Finally, the Airy 2 process occurs in equilibrium statis- 
tical mechanics too. It describes the fluctuation of the border of a facet in the 
3d-Ising corner model [19,38], but the same result is expected to hold in much 
more generality for short-range interactions crystals at low temperature [18], 
see also [33] for related studies on the macroscopic shapes. 

We have seen that the GUE Tracy- Widom distribution and the Airy 2 pro- 
cess appears in several models. How stable are they under change of initial 
conditions or intrinsic symmetries of the model? While the scaling exponents 
remain universal, the precise distribution depends on initial conditions or sym- 
metries of the model. However, the results obtained so far indicate that the 
limit distribution is still universal but one has to divide the universality class in 
a few subclasses. 

In classical random matrix theory, besides the GUE, there are the Gaus- 
sian Orthogonal Ensemble (GOE) and the Gaussian Symplectic Ensemble 
(GSE) [15]. They were introduced as an approximation to the real Hamilto- 
nian of heavy nuclei and the different ensembles reflect the different intrinsic 
symmetries of the system under consideration (e.g., the time- reversal symme- 
try, which is broken by external magnetic fields). With the same scaling as for 
GUE, Tracy and Widom obtained the limit distribution of the largest eigen- 
values for GOE and GSE too [51,53]. This is a clear example of the stronger 
stability of the scaling exponent (in this case the fluctuation exponent) v.s. the 
distribution functions. In the PNG model, the GOE Tracy- Widom distribution, 
F\ 1 occurs when the growth is on a flat substrate [5,40]. The connection be- 
tween the PNG model and the GOE ensemble actually goes beyond the largest 
eigenvalue [16]. 

Since 2002, the question of the analogue of the Airy 2 process in this setting 
was open and, for some aspects, like in random matrix, is still open. In 2005 
Sasamoto discovered a new process in the context of the TASEP [44], see [10] 
for a complete derivation. This process is now called the Airyi process and, 
although has a similar mathematical structure as the Airy 2 process, a lot of the 
probabilistic interpretation that was present for the latter is lost (or, maybe, not 
yet understood). The first generalization of the setting considered in [44] has 
been carried out in [9] . The evidence of universality of the Airyi process is not 
yet as large as for the Airy 2 process, but it is expected to occur in point-to-line 
directed last passage percolation, in the polynuclear growth model on a flat (not 
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necessarily horizontal) substrate, and possibly also in GOE Dyson's Brownian 
Motion. 

In this paper we first define the Airyi and Airy2 processes and give their 
known properties. Then we introduce the continuous time TASEP and explain 
under which initial conditions the two Airy processes occurs. For the interested 
reader around some of the methods used in relation to the Airy2 process as well 
as related models, we refer to a the lecture notes [31,49], and the surveys [12, 
17,24,35,48]. 



2 The Airyi and Airy2 processes 

In this section we define the two Airy processes and give the known properties. 

Definition 2.1 (The Airyi process). The Airyi process A\ is the process with 
m-point joint distributions at U\ < u 2 < ■ ■ ■ < u m given by the Fredholm deter- 
minant 

m 

P ( f] Ml(«k) < s k}) = det(l - XsK 1 Xs)L^{u 1 ,...,u m }xR) (2.1) 

fe=l 

where Xs(uk,x) — l(x > s^) and the kernel K\ is a follows. Let 

A\(x, y) — Ai(x + y) and H\ = —A with A the one- dimensional Laplacian. 
Then the K\ is defined by 

Ki(u, s; u', a') = -(e~^ u '-^ Hl )(s, s')l(u < u') + (e uHl A ie - u ' Hl ){s, s'). (2.2) 

Definition 2.2 (The Airy2 process). The Airy2 process Ai is the process with 
m-point joint distributions at ui < u 2 < . . . < u m given by the Fredholm deter- 
minant 

m 

P( f) {A 2 (u k ) < = det(l - XsK 2X s)m{u u ...,u m }xR) (2.3) 
fe=i 

where Xs{uk,x) = l(x > Sk) and K 2 is the extended Airy kernel. Let H 2 be the 
Airy operator 

< 2 - 4 > 

and A 2 the (one-time) Airy kernel with entries 

A 2 (x,y)= [ dAAi(x + A)Ai(y + A). (2.5) 

Then the define the kernel 
K 2 (u,s;u',s') = -(e-("'-") ff2 )( S , S ')l(« <u') + {e uH iA 2 e- u ' H *)(s,s'). (2.6) 

Explicit expressions 

The explicit formulas of the K\ and K 2 kernels are the following. For the Airyi 
process, as shown in Appendix A of [10], one has 



+Ai(s + s + (u - uf) cxp ([v! - u)(s + s') + |(u' - uf^j 



(2.7) 
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with Ai the Airy function [1], while the kernel for the Airy2 process writes [29,42] 

f L dAe("'~") A Ai(x + A)Ai(y + A), u > u', 
K 2 (u,s;u',s') = \ R+ , , „ " (2.8) 

\ -/„ dAe( u - tl ) A Ai(a; + A)Ai( 2 / + A), u < u'. 

One-point distributions 

The Airyi and Airy2 processes are stationary and their one-point distribution 
functions are given in tems of the GOE and GUE Tracy- Widom distributions 
of random matrices [50,51]. Namely, as shown in [20], 

P(A(0)<s)=Fi(2s), (2.9) 

and, for the Airy2 process [42], 

W(M0) < s) = F 2 (s). (2.10) 

Spatial correlations 

Some information on the spatial correlation have been already determined for 
the Airy2 process. Locally it behaves like a diffusion and on a long distance the 
spatial correlations have a slow polynomial decay. More precisely, define the 
function g by 

Var(A 2 (u)-A 2 (0))=g(u). (2.11) 

g grows linearly for small u and that the Airy process has long range correla- 
tions [42]: 

- / 2u + °( u2 ) for M Small > (0 10\ 

9{ > ~ \ g(oo) - 2u- 2 + 0(u- 4 ) for |u| large. 1 > 

with g(oo) — 1.6264 . . .. The coefficient 2 of the correlation's decay is determined 
in [2,54]. For the Airy 2 process a set of PDE's [2] and ODE's [52] have been 
obtained. 

For the Airyi process such analysis is still missing, but it is expected that 
the global behavior of g is the same, as indicated also from some preliminaly 
numerical computations [44]. 

On the Fredholm determinants 

The joint distributions for the Airy processes are given in terms of Fredholm 
determinants. Are they really well defined? One way of considering the Fred- 
holm determinants (|2.ip and (|2.3|) is simply via their Fredholm series expansion, 
namely 

det(l - XsK X s)lz « m }xR) 

(2.13) 

(—1)™ "\ r r 
= ^ ~ X! / & x i---j da;„det(K(u i(c ,a;fc;ui i ,a;())i<fe,i = 



Kn- 



One can see that (|2.13j) is absolutely summable/integrable for the kernel K = K\ 
and K = K^. However, it is sometimes useful to consider the expressions 
(|2.ip and (|2.3p as Fredholm determinant of an operator on the Hilbert space 
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rate 1 , , rate 1 




Figure 1: Dynamics of the TASEP. Particles jump with rate one on the right 
site, but under the constraint that the site is empty. 



H = L 2 ({ui, . . . , u m } x K). The Fredholm determinant is well defined for trace 
class operators. 

For (|2.3|) there is no problem, since is trace class on Ji [50]. This is not 
the case for K\. In fact, due to contribution of the diffusion part of the kernel, 
K\ is not even Hilbert-Schmidt on H. However, as shown in Appendix A of [9], 
there exists a conjugate operatorof K\ that is trace class on TL. With conjugate 
operator we mean an operator K\ that leads to the same correlation functions 
and Fredholm series expansion, i.e., such that 

det(Ki(u lk , x k ; u H , xi))i< k ,i< n = det(Ki(u lk , x k ; u it , Xi))i<k,i< n (2-14) 

holds. Equivalently, one can employ a weighted Hilbert space where K\ is trace 
class. 

3 The TASEP 

The totally asymmetric simple exclusion process (TASEP) is an interacting 
stochastic particle system. It consists on particles on Z with the exclusion con- 
straint that at any given time t, every site j G Z can be occupied at most 
by one particle. Thus a configuration of the TASEP can be described by 
r\ = {r)j,j G Z|?7j G {0, 1}}. rjj is called the occupation variable of site j, which 
is defined by rjj — 1 if site j is occupied and r\j — if site j is empty. 

The dynamics of the TASEP is defined as follows. Particles jumps on the 
neighboring right site with rate 1 provided that the site is empty, see Figure [TJ 
This means that jumps are independent of each other and are performed after an 
exponential waiting time with mean 1, which starts from the time instant when 
the right neighbor site is empty. More precisely, let /: SI — > M be a function 
depending only on a finite number of r/j 's. Then the backward generator of the 
TASEP is given by 

Lf(v) =E^'( 1 -^+ 1 )W' i+1 )-/W)- (3- 1 ) 

Here rf'^ 1 denotes the configuration rj with the occupations at sites j and 
j + 1 interchanged. The semigroup e Lt is well-defined as acting on bounded and 
continuous functions on fi. e Lt is the transition probability of the TASEP [34]. 

Observables 

Since the dynamics does not interchange particles, they can be labelled by an 
index k G /, I C Z. We denote by Xk(t) the positions of the particle k at time 
t and use the right-left ordering, i.e., Xk+i(t) < x k {t), for all t. There are two 



5 



observables (quantity of interest) which are closely related: the position of given 
particles and the integrated current at fixed locations. 

The TASEP integrated current at position x and time t, J(x, t), is the num- 
ber of particles which jumped from site x to site x + 1 during the time interval 
[Q,t\. Let us label by 1 the right-most particle starting at position x\(0) < x. 
Then J(x,t) and x s (t) are related by 

¥(J(x, t)>s) = ¥(x s (t) > x + 1). (3.2) 

Below we consider as observable a given subset of particles. For some finite 
/ C Z, we explain in which scaling limit the joint distributions of 

{x k (t),kel} (3.3) 

are governed by the Airy processes. 



Step initial conditions and Airy 2 process 

The Airy2 process appears, is called step-initial condition. The initial configu- 
ration is deterministic defined by 

x k (t = 0) = -k, ke 1+. (3.4) 

First consider the macroscopic behavior, around which one analyzes the fluctu- 
ations. The macroscopic limit density u(£) is given by 

u(0 = A lim t" 1 E(#(fc : x k (t) > &)). (3.5) 

It was proven by Rost [43] that the macroscopic density has a linearly decreasing 
region, namely 

r i, €<-i, 

u(0 = { 1 - l)/2, 1 < C < 1, (3.6) 
I 0, £>1. 

The particles which build the bulk of the linear region are the one with parti- 
cles number at, a G (0, 1). First it was shown, by means of a growth model [27], 
that the fluctuations of their position is F 2 -distributed. Then Johansson in [29] 
proves a functional limit theorem in a discrete-time setting. Its continuous-time 
version writes 

Urn '[■/^■"iW^y-'C/"^ . Mu) . (3 . 7) 

t^oo — (t/2) 1 / ci 

Alternating initial conditions and Airyi process 

The situation in which the Airyi process was first discovered is the alternating 
initial condition. It is the deterministic initial configuration given by 

x k (t = 0) = -2k, keZ. (3.8) 

The macroscopic density is simply u(£) = 1/2 and particles moves with 
average speed 1/4. For alternating initial condition, it was expected by univer- 
sality that the one-point distribution was Fi-distributed. It was the case for a 
continuous version of a corresponding growth model [4, 40] . 
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The extension to joint distribution remained for a while unsolved. The 
solution came from a new approach. Starting from a determinantal formula by 
Schiitz [45] of the joint distribution of particle positions, Sasamoto [44] found 
a clever way of rewriting in term of a signed determinantal point process on a 
larger set of variables. This is explained in details in [10], where it is also proven 
the pointwise convergence of kernel to K\ . The finite-dimensional convergence 
to the Airyi process has been obtained for a discrete-time version of the TASEP 
in [9]. Its continuous-time analogue writes 

t ^ Z?73 = Mu)- (3.9) 



Stationary TASEP 

Besides deterministic initial conditions, one can also consider the stationary 
TASEP. The translation invariant measures for the stationary TASEP are the 
Bcrnoulli-p measures, p £ [0,1], i.e., P(% = 1) = p and the random variables 
?7fc's are independent. Between other quantities, an explicit expression for the 
one-point distribution has been obtained in continuous time and discrete time 
TASEP with parallel update [4,22,41]. 



Universality and generalizations 

By universality it is expected that the Airyi process is the large time limit pro- 
cess for a larger class of models. First of all, the Airyi process should appear 
for all deterministic and periodic initial conditions for the TASEP, and its dis- 
crete time versions. It is also expected that the Airy processes occur in the 
partially asymmetric simple exclusion process (PASEP), as soon as the drift is 
non-zero. As it will be explained below, the TASEP can be reinterpreted as a 
growth model belonging to the KPZ universality class. For model in this class of 
universality, we expect that the Airy processes appears too. More precisely, the 
Airyi process for growth on a flat substrate and the Airy2 process for growth 
leading a curved shape due to initial conditions. Some of these growth models 
already studied are directly linked with directed last passage percolation. Thus, 
we expect the Airyi, resp. Airy2, process to describe the process associated to 
point-to-linc, resp. point-to-point, directed last passage percolation. 

As discussed briefly in the Introduction, the Airy 2 process has been obtained 
in several models. For the Airyi process only in a few cases the analysis is 
available so far. The first generalizations has been made in [9], where the TASEP 
is analyzed for a larger set of initial conditions and in the discrete-time TASEP 
with sequential update. This update rule is as follows. At each time step 
particles are processed sequentially from right to left, i.e., starting from right 
to left, if the site on the right of a particle is empty, then this particle jumps 
there with probability p. This update rule allows to shift blocks of particles to 
the right during one time-step. 

More precisely, in [9] it is considered the deterministic and periodic initial 
configuration where particles start from dZ for any d > 2 fixed. Then, it is 
proven that the particle positions converges, properly rescaled, to the Airyi 
process. More precisely, 

,. x n{ ^ t) (t) - put 2 / 3 

t 1 ™ Z^fs = Mu), (3.10) 
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where the constants /i and k are given by 



wss (2(l-p) P )V»(d(d-l))V' ^ 2 

a — p a — 1 



and the index of the particle n(u, t) by 
n(u, t) 



P(d- l) f _ M^2/3 
d(d — p) d 



(3.12) 



The convergence is in the sense of finite-dimensional distributions. 

In discrete time, a second natural update rule is the parallel update. It 
consists in first checking for all particles if they can jump (i.e., if their right- 
neighbor is empty) and then, simultaneously and independently, these particles 
jump to the right with probability p. Some progress in the parallel update has 
been recently made in [39], but there are not yet results concerning the limit 
processes. Other update rules can also be considered, see the review [46]. 

A natural extension is to consider the partially asymmetric simple exclusion 
process (PASEP), but very few is known. The scaling exponent have been 
determined also on a rigorous level [7] , but most of the methods applied for the 
TASEP do not easily extend to the PASEP. Also in this model the Airyi and 
Airy2 processes are expected to occur as soon as the drift is non-zero. 



A geometric representation 

In this final part we present a geometric representation of the TASEP, from 
which the reason why the Airyi process should occur in the point-to-line last 
passage directed percolation is apparent. 

Growth model. Given a configuration of particles, one associates a height func- 
tion h. Set the height at a given position, e.g., h(Q) = 0. Then, the height 
function at x g Z is obtained as follows. The height differences are given by the 
occupation variables of the TASEP, more precisely, h(x + 1) — h(x) = 1 — 2-q x . 
Thus, starting from x = 0, one can define the height h{x), ifZ. The extension 
to x G R is just by linear interpolation of the heights on Z. An example is shown 
in Figure O 




x 




Figure 2: The interface (height line) associated with a TASEP configuration. 



The dynamics of the TASEP is reflected in the height function picture by a 
stochastic growth in the vertical direction. In fact, when a particle jumps from 
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x to x + 1, in the height function a valley \/ becomes a mountain /\ and the 
height at x + 1 is increases by two, see Figure [5] In this way, the TASEP can 
be equivalcntly seen as a stochastic growth model of an interface. 

The two types of initial conditions for the TASEP discussed above corre- 
sponds to the so-called corner and flat growth. More precisely, step initial con- 
ditions become growth starting from a corner-like initial surface, while periodic 
initial conditions become growth from a flat (also tilted) initial configuration, 
see Figure [31 With flat is it meant that the surface is flat after a coarse graining 
of order one. 




Figure 3: Interfaces associated with a step initial configuration (left) and two 
periodic initial configurations (right). 



Corner growth leads then to a macroscopic shape which is curved and the 
large time limit process describing the surface is the Airy2 process, while growth 
on a flat substrate keeps a macroscopically flat interface with fluctuations de- 
scribed by the Airyi process. 

For discrete TASEP with parallel update, by setting space and time units 
to p and letting p — > 0, one obtains the (continuous time) polynuclear growth 
(PNG) model [40,42], a stochastic growth model in the KPZ class. 

Directed last passage percolation. The TASEP can be mapped to last passage 
percolation on Z 2 with i.i.d. exponentially distributed random variables u>(m, n), 
m,n G Z. The precise connection is that co(m, n) is the waiting time of particle 
number n to jump from position to — n — 1 to m — n. 

To explain it better, consider the step-initial condition, which, as we will 
see, corresponds to point-to-point directed last passage percolation. The last 
passage time G(to, n) is given by 

G(m,n)= max £(ir), £(ir) = V w(i,j) (3.13) 

7r:(l,l)->(m,n) f-f 

where n : (1,1) — » (to, n) are up-right paths (concatenation of (0,1) or (1,0) 
steps), going from (1, 1) to (m, n). If there are no up-right paths from (1, 1) to 
(to, 71) , we set G(m,n) = 0. Then define the domain 

B t = {(m,n) e Z 2 \G(m,n) < t}. (3.14) 

The border of this domain is directly connected with the TASEP as follows. 

To the step-initial configuration 2^(0) = — fc, k — 1,2,..., we associate 
the height line of Figure [3] and after a clockwise rotation of 45 degrees, one 
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obtain the line of Figured] (left). This line is the border of Bo, denoted by 
dBo, which passes through the coordinates (xk(0) + k, k), k = 1,2,.... Then, 




Figure 4: TASEP with step-initial condition and directed percolation. Left is 
the £ = configuration and right is the one at time £ = w(l, 1). 



at time £ = a; (1,1), dBt becomes the one of Figure H] (right). This happens 
exactly when the particle number 1 jumps from position —1 to position 0. The 
correspondence actually extends for all £ > and using the identity G(m, n) = 
w(m, n) + max{G(m — 1, n), G(m, n — 1)} one can see that G(m, n) is the time 
needed for particle number n to reach site m — n. 

Finally, periodic initial conditions become the point-to-line problem in last 
passage percolation. More precisely, the tt in (|3.13[) are all paths ending at (m, n) 
and starting from a given line, e.g., the line is {(it, —u), u G Z} for Xfc(0) = — 2fc, 
k e Z. 
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